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Abstract: We show that the entropy created by Ohmic dissipation inside an accreting 
charged black hole may exceed the Bekenstein-Hawking entropy by a large factor. If the black 
hole subsequently evaporates, radiating only the Bekenstein-Hawking entropy, then the black 
hole appears to destroy entropy, violating the second law of thermodynamics. A companion 
paper discusses the implications of this startling result. Bousso's covariant entropy bound is 
not violated. 
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1. Introduction 

The purpose of this paper is to show that classical processes of dissipation can generate huge 
quantities of entropy inside the horizon of a black hole, many orders of magnitude more than 
the Bekenstein-Hawking Q entropy. The specific black hole model presented is intended to 
be semi-realistic, albeit over-simplified, with parameters appropriate to a real supermassive 
black hole. We take charge as a surrogate for angular momentum, and electrical conductivity 
as a surrogate for angular momentum transport. To see how much entropy might be created, 
we treat the electrical conductivity as an adjustable free parameter. 

If a black hole creates many times the Bekenstein-Hawking entropy and subsequently 
evaporates, radiating only the Bekenstein-Hawking entropy and leaving no remnant, then 
entropy is destroyed, violating the second law of thermodynamics. The implications of this 
startling result are discussed in a companion paper Q. 

Throughout this paper we treat entropy in a purely classical fashion. In particular, we 
assume that locality holds inside the black hole. Locality, the quantum field theory proposition 
that operators commute at spacelike-separated points, is the assumption that normally makes 
it legitimate to add entropy over spacelike surfaces. Since the spacetime curvature inside a 
supermassive black hole is well below Planck, except near the singularity, one might expect 
classical physics to apply. 

It is widely thought that in order to preserve unitarity of black hole evaporation, locality 
must break down over spacelike surfaces connecting the inside and outside of a black hole || . 
In the companion paper ||] we argue that the calculations of the present paper point to a 
profligate breakdown of locality inside black holes. 

We work in Planck units, ks = c = G = h = 1. 
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2. Model of entropy production inside a black hole 

Real supermassive black holes acquire most of their mass not during a single collapse event, 
but rather by gradual accretion of gas over millions or billions of years. We model this 
gradual growth by the general relativistic, self-similar, accreting, spherical, charged black hole 
solutions described by [|], g], to which the reader is referred for more detail. In these models 
the black hole accretes a charged "baryonic" plasma of relativistic matter, with constant 
proper pressure-to-density ratio p/p = w, at a constant rate, so that the mass of the black 
hole increases linearly with time. 

Real supermassive black holes probably rotate, but have tiny electric charge. However, 
the interior structure of a spherical charged black hole resembles that of a rotating black hole in 
that the negative pressure (tension) of a radial electric field produces an effective gravitational 
repulsion analogous to the centrifugal repulsion inside a rotating black hole. Thus we follow 
the common practice of taking charge as a surrogate for spin. In the self-similar solutions, 
the charge of the black hole is produced self-consistently by the accumulation of the charge 
of the accreted plasma. 

Similarly, we take electrical conduction as a surrogate for the dissipative transport of 
angular momentum. We do not use a realistic electrical conductivity, but rather treat it as 
a phenomenological adjustable quantity. In diffusive electrical conduction, the electric field 
E = Q/r 2 gives rise to an electric current j = aE. If the conductivity a is taken to be a 
function only of the plasma density p, then the condition of self-similarity forces Q 

a = kp 1 I 2 /{Att) 1 / 2 (2.1) 

where the dimensionless conductivity coefficient k is a phenomenological constant. As dis- 
cussed by Q, this phenomenological conductivity is greatly suppressed compared to any 
realistic conductivity (except perhaps at densities approaching the Planck density). How- 
ever, angular momentum transport is intrinsically a much weaker process than electrical 
conduction, so it is not unreasonable to consider a greatly suppressed conductivity. 

Since information can propagate only inwards inside a black hole (at least classically) , it is 
natural to impose boundary conditions outside the black hole. We assume that the boundary 
conditions of the accreting black hole are established at a sonic point outside the horizon, 
where the infalling plasma accelerates from subsonic to supersonic. We assume that the 
acceleration through the sonic point is finite and differentiable, which imposes two boundary 
conditions. The accretion in real black holes is likely to be much more complicated, but this 
assumption is the simplest physically reasonable one. 

We define the charge Q, and mass M, of the black hole at any instant to be those that 
would be measured by a distant observer if there were no charge or mass outside the sonic 
point, 

O 2 

Q. = Q and M. = M + -j— at the sonic point (2.2) 
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where r is the circumferential radius, and Q and M denote the interior charge and mass, all 
gauge-invariant scalar quantities. If the black hole ceases accreting abruptly at some time, 
then Q m and M, are the actual charge and mass of the black hole at that time. 

Given the assumption that the sonic point is regular, the dimensionless free parameters 
of the solutions are: (1) the mass accretion rate M»; (2) the charge-to-mass ratio Q,/M, of 
the black hole; (3) the equation of state parameter w; and (4) the conductivity coefficient k. 

The black hole mass increases linearly with time, M, oc t, and the mass accretion rate 
M, is 

M. = dM./dt = M./t , (2.3) 

where t = Td = {r^ t d ) r=Ta is the time measured by clocks attached to neutral dust (d) that 
free-falls radially through the sonic point r = r s from zero velocity at infinity, and which 
therefore records the proper time at rest at infinity, and £5 is the time component of the 
homothetic vector £ k in the dust frame |], Appendix E]. 

The density p and temperature T of a relativistic fluid in thermodynamic equilibrium 
are related by p = (7r 2 g/30)T 4 , where g = gs + |sf is the effective number of relativistic 
particle species, with gs and gp being the number of bosonic and fermionic species. If the 
expected increase in g with temperature T is modeled (so as not to spoil self-similarity) as a 
weak power law g/g p = T £ , with g p the effective number of relativistic species at the Planck 
temperature, then the relation between density p and temperature T is 

p = (vr V30)T( 1+ -)/- , (2.4) 

with equation of state parameter w = 1/(3 + e) slightly less than the standard relativistic 
value w = 1/3. We fix g p by setting the number of relativistic particles species to g = 5.5 at 
T = lOMeV, corresponding to a plasma of relativistic photons, electrons, and positrons. 

The entropy S of a proper Lagrangian volume element V of an ideal relativistic fluid with 
zero chemical potential is S = [(p + p)/T]V. The proper velocity of the baryonic fluid through 
the sonic point equals the ratio £ r /£* of the radial and time components of the homothetic 
vector in the plama frame Q . Thus the entropy S accreted through the sonic point per unit 
proper time of the fluid is dS/dr = [(1 + w)p/T]4irr' 2 (£ i r /£*). The sonic radius r s of the black 
hole increases as d\nr s /dr = l/(r^ t ) r=ra [Q]. The Bekenstein-Hawking entropy of the black 
hole is 5bh = tit^ where is the horizon radius, so cLSbh/ 'd\nr s = Thus the entropy 

S accreted per unit increase of the Bekenstein-Hawking entropy is 

dS _ _L_ (1 + w)piirr 3 e 
dS BU ~ 2vrr2 T r=Ta ' ( ' > 

Inside the sonic point, dissipation increases the entropy. The energy-momentum tensor 
is the sum of baryonic and electromagnetic parts, T™ n and T™ n , and the evolution of baryon 
entropy is determined by the time component of the equation of covariant conservation of 
energy-momentum in the rest frame of the baryons: 

D m T^ = -B m Tf . (2.6) 
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In the self-similar model being considered, the energy conservation equation ( |2.6| ) can be 
written [01 



dp riln(r 3 Q Q 2 



(2.7) 



which can be recognized as an expression of the first law of thermodynamics dpV+pdV = TdS 
with proper volume V oc r 3 £ r . The right hand side of equation (2/7) is the Ohmic dissipation 
jE = gE 2 . Equation (|2.7|) can be re-expressed as 



din 5 
dr 



aQ 2 



(1 +w)pr 4 



{21 



with S oc /0 1 /(l+») r 3^ (x ( /0 /T)r 3 e r . 

Since other physics presumably takes over near the Planck scale, we truncate the pro- 
duction of entropy at some arbitrary maximum density put ( "rho splat" ) . Integrating equa- 
tion ( |2.8| ) from the sonic point to the splat point yields the ratio of the entropies at the sonic 
and splat points. Multiplying the accreted entropy, equation (gj^), by this ratio yields the 
rate of increase of the entropy of the black hole, truncated at the splat point, per unit increase 
of its Bekenstein-Hawking entropy 



dS 
gLSbh 



1 (1 +w)pATrr 3 C 



2nr 2 



(2- 



p=p# 



The entropy created becomes large when the conductivity coefficient lies in the range 
k ps 1.3 to 1000. Over this range the rate dS/dS-sn of increase of entropy, equation (|2.9D , is 
almost independent of the black hole mass M, , 



dS 
dSBu 



const 



2(l-Q2 /M 2)i/2 (1 + w)/9# 
M.[1 + (1-Q2/M. 2 )V 2 ]V # T # 



(2.10) 



in which the empirical fit on the right hand side is accurate to a factor of two over the 
range k ~ 10 to 1000 (for k < 10 to ~ 1.3, the fit increasingly overestimates dS/dS^n), and 
M, > 3M Q , M. < 10~ 4 , Q./M. ps 10" 12 to 0.99999, w m 0.1 to 0.55, and p # not too small. 

Bousso 0] has proposed the covariant entropy bound, which states that the entropy 
passed through a converging lightsheet cannot exceed S cov = area/4 of its boundary. In the 
models under consideration, the rate at which entropy passes through an ingoing or outgoing 
spherical lightsheet per unit decrease in S a 



nr 2 is 



dS 



dS c 



dS 



dS BU r 2 ei/3T7l 



(2.11) 



in which {/3,7} are the time and radial components of the proper covariant radial 4-gradient 
in the notation of [0J, and the =F sign is — for ingoing, + for outgoing lightsheets. A sufficient 
condition for the covariant entropy bound to be satisfied is \dS/dS cov \ < 1. 
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Figure 1: (Color online) Internal structure of a black hole with mass M. = 4x 10 6 M Q , accretion 
rate M. = 1CP 16 , charge to mass Q,/M, = 10~ 5 , equation of state w — 0.32, and conductivity 
coefficient k = 1.24. The quantities plotted function of radius r: the density p of the baryonic 

plasma, the energy density p e (short dashed line) of the static electric field, the absolute value of the 
Weyl curvature scalar C = Airp/S + Q 2 / (2r 4 ) — M/r 3 , the rate dS/dSBa of increase of the black hole 



entropy with Bekenstein-Hawking entropy, equation (2J3), and the rate \dS/dS cov \ at which entropy 
passes through ingoing (dark) and outgoing (light) spherical lightshects per unit decrease in their 



area/4, equation (2.11). Vertical dashed lines mark the horizon, and where the Weyl curvature \C\ 
exceeds 1 Planck unit. Arrows, such as that above dS/dSBH, show how the curves shift a factor of ten 
into the past and the future. The rate dS/dSBn is almost independent of the black hole mass M,, at 
fixed splat density p#. 



3. Example of interior entropy exceeding Bekenstein-Hawking 

A black hole of mass 4 x 10 6 M « 4 x 10 44 Planck units (the mass of the supermassive 
black hole at the center of the Milky Way ||, ||]) accreting over the age of the Universe 
10 10 yr ~ 6 x 10 60 Planck units has an accretion rate of M, ~ 1CP 16 . Figure |l] shows 
the interior structure of a black hole with that mass and accretion rate, charge-to-mass 
Q 9 /M m = 10 -5 , equation of state w = 0.32, and conductivity coefficient k = 1.24. 

To produce lots of entropy, the baryonic plasma must fall to a central singularity, and we 
choose the conductivity k = 1.24 to be at (within numerical accuracy) the critical conductivity 
[U for this to occur. Below the critical conductivity the plasma generally does not fall to a 
singularity, but rather drops through the Cauchy horizon. The latter solutions are subject to 
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the mass inflation instability [|i"o"| , ||, a fascinating regime not considered in this paper. 

Solutions at the critical conductivity exhibit § the periodic self-similar behavior first 
discovered by [O]. The ringing of the curves shown in Figure || is a manifestation of this, not 
a numerical error. 

The electric charge advected by the plasma inwards across the sonic point is, thanks 
to the "high" conductivity, almost canceled by an outward current. As a consequence, the 
charge-to-energy of the accreted plasma, here « 400 at the sonic point, is substantially larger 
than the charge-to-mass of the black hole. We have deliberately chosen a small charge-to- 
mass ratio for the black hole, Q,/M, = 10 -5 , so that the Lorentz repulsion of the plasma by 
the black hole is subdominant, and the trajectories of parcels of plasma outside the black hole 
are not greatly different from Schwarzschild geodesies. Thus for example the sonic point is at 
a radius of 3.06 geometric units (c = G = M. = 1), close to that expected for a neutral tracer 
relativistic fluid that free-falls from zero velocity at infinity. The horizon is at 2.00 geometric 
units, like Schwarzschild. Figure |l] shows the solution out to 2,000 geometric units. 

At the sonic point, the plasma temperature is ~ 4 x 10 5 K. Inside the horizon, the electric 
field increases, and Ohmic dissipation starts to heat the plasma, increasing its temperature 
and entropy. When the plasma energy has become comparable to the electric energy, then 
the plasma goes into a power law regime where the plasma and electric energies increase in 
proportion to each other, kept in lockstep by the conductivity. 

The entropy hits the Bekenstein-Hawking milestone, dS/dSsu = 1, when the temperature 
is ~ 3 x 10~ 16 Planck units, or 3TeV, and the curvature radius is |Cj -1 / 2 ~ 10 30 Planck 
lengths, or 0.01 mm. This temperature and curvature are almost independent of the mass M, 
of the black hole, equation ( 2.1 0| ) . 

If the plasma's dissipative trajectory is followed to the Planck scale, \C\ = 1, then the 
rate of increase of entropy relative to Bekenstein-Hawking is dS/dSsu ~ 10 10 , again almost 
independent of the mass M, of the black hole. If the entropy is assumed to accumulate 
additively inside the black hole then the cumulative entropy can evidently exceed Bekenstein- 
Hawking by a large factor. 

Figure |l] shows that \dS/dS cov \ < 1, equation ( 2.11 ), at all sub-Planck scales. Thus 
although the cumulative entropy may exceed Bekenstein-Hawking, Bousso's covariant entropy 
bound is satisfied by the black hole. 



4. Conclusion 



We have shown that the dissipation of the free energy of the electric field inside a charged 
black hole can potentially create many times more than the Bekenstein-Hawking entropy. If 
the black hole subsequently evaporates, radiating only the Bekenstein-Hawking entropy and 
leaving no remnant, then entropy is destroyed. 

This startling conclusion is premised on the assumption that entropy created inside a 
black hole accumulates additively on spacelike slices, which in turn derives from the assump- 
tion that the Hilbert space of states is multiplicative over spacelike-separated regions, as 
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postulated by locality. This is essentially the same reasoning that originally led Hawking 
[p^ l to conclude that black hole evaporation is non-unitary, that black holes must destroy 
information. 

It is widely thought that unitarity should be considered a higher principle than locality. 
To ensure that black hole evaporation is unitary, locality between the inside and outside of 
a black hole must break down. In a companion paper j|] we argue that the gross violation 
of the second law found in the present paper points to a wholesale breakdown of locality 
inside black holes, and provides a compelling argument in favor of the conjecture of "observer 
complementarity" . 

The black hole respects Bousso's covariant entropy bound J5J, as it should given the 
theorem of [O]. 
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